function out = L1Energy(u, g, lambda)
y = NeuGradient(u);
s = 0;
s = sum(sum(abs(y))) + lambda * sum(sum(abs(u - g)));
out = s;
end
